Adaptive synchronization in delay-coupled networks of Stuart-Landau oscillators 
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We consider networks of delay-coupled Stuart-Landau oscillators. In these systems, the coupling 
phase has been found to be a crucial control parameter. By proper choice of this parameter one can 
switch between different synchronous oscillatory states of the network. Applying the speed-gradient 
method, we derive an adaptive algorithm for an automatic adjustment of the coupling phase such 
that a desired state can be selected from an otherwise multistable regime. We propose goal functions 
based on both the difference of the oscillators and a generalized order parameter and demonstrate 
that the speed-gradient method allows one to find appropriate coupling phases with which different 
states of synchronization, e.g., in-phase oscillation, splay or various cluster states, can be selected. 
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I. INTRODUCTION 



The ability to control nonlinear dynamical systems has 
brought up a wide interdisciplinary area of research that 
has evolved rapidly in the last decades [l[. In particular, 
noninvasive control schemes based on time-delayed feed- 
back [2|t4| have been studied and applied to various sys- 
tems ranging from biological and chemical applications to 
physics and engineering in both theoretical and experi- 
mental works (5l-[ia|. Here we propose to use adaptive 
control schemes based on optimizations of cost or goal 
functions |13l4l5l | to find appropriate control parameters. 
Besides isolated systems, control of dynamics in spatio- 
temporal systems and on networks has recently gained 
much interest |16l - [2(j | . The existence and control of clus- 
ter states was studied by Choc et al. [2l|, [g^] in networks 
of Stuart-Landau oscillators. This Stuart-Landau sys- 
tem arises naturally as a generic expansion near a Hopf 
bifurcation and is therefore often used as a paradigm for 
oscillators. The complex coupling constant that arises 
from the complex state variables in networks of Stuart- 
Landau oscillators consists of an amplitude and a phase. 
Similar coupling phases arise naturally in systems with 
all-optical coupling [a, [2^. Such phase-dependent cou- 
plings have also been shown to be important in over- 
coming the odd- number limitation of time-delay feedback 
control [2J, |25| and in anticipating chaos synchronization 
|26| . Furthermore, it was shown in [2l|,[23l that the value 
of the coupling phase is a crucial control parameter in 
these systems, and by adjusting this phase one can delib- 
erately switch between different synchronous oscillatory 
states of the network. In order to find an appropriate 
value of the coupling phase one could solve a nonlinear 
equation that involves the system parameters. However, 
in practice the exact values of the system parameters are 
unknown, and analytical conditions can be derived only 



for special values of the complex phase. An efficient way 
to avoid these limitations and find optimal values of the 
coupling phase is the use of adaptive control. 

In this paper, we present an adaptive synchronization 
algorithm for delay-coupled networks of Stuart-Landau 
oscillators. To find an adequate coupling phase we apply 
the speed-gradient method [15| . which was used previ- 
ously in various nonlinear control problems, yet not for 
the control of dynamics in delay-coupled networks. By 
taking an appropriate goal function we derive an equa- 
tion for the automatic adjustment of the coupling phase 
such that the goal function is minimized. At the same 
time the coupling phase converges to the theoretically 
predicted value. Our goal function is based on the Ku- 
ramoto order parameter and is able to distinguish the 
different states of synchrony in the Stuart-Landau net- 
works irrespectively of the numbering of the nodes. 

This paper is organized as follows. After this introduc- 
tion, we describe the model system in Sec. [Ill Section Hill 
introduces the speed-gradient method and its application 
using the coupling phase in networks of Stuart-Landau 
oscillators. We present the main results for the control of 
in-phase synchronization in Sec. IIVI and for cluster and 
splay states in Sec. |Vl Finally, Sec. IVII contains some 
conlclusions. 



II. MODEL EQUATION 

Consider a network of N delay-coupled oscillators 

N 



with Zj ~ rjc'-'fi £ C, j ^ 1, . . . ,N. The coupling matrix 
A = {aij}fj^i determines the topology of the network. 
The local dynamics of each element is given by the normal 
form of a supercritical Hopf bifurcation, also known as 
Stuart-Landau oscillator. 
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f{z,) = [X + ioj-{l+ij)\zj\^]zj 



(2) 



with real constants A, a; 7^ 0, and 7. In Eq. ([T]), t is the 
delay time. K and j3 denote the amplitude and phase of 
the complex coupling constant, respectively. Such kinds 
of networks are used in different areas of nonlinear dy- 
namics, e.g., to describe neural activities |27| . 

Synchronous in-phase, cluster, and splay states are 
possible solutions of Eqs. ([T]) and ^. They exhibit 
a common amplitude tj = ro,m and phases given by 
ipj — Qmt + jAipm with a phase shift Aipm = 2'Km/N 
and collective frequency flm- The integer m determines 
the specific state: in-phase oscillations correspond to 
771 = 0, while splay and cluster states correspond to 
m = 1,...,A^ — 1. The cluster number d, which de- 
termines how many clusters of oscillators exist, is given 
by the least common multiple of 777 and TV divided by 771, 
and d = N (e.g., 777 = 1), corresponds to a splay state. 

The stability of synchronized oscillations in networks 
can be determined numerically, for instance, by the mas- 
ter stability function [28| . This formalism allows a sep- 
aration of the local dynamics of the individual nodes 
from the network topology. In the case of the Stuart- 
Landau oscillators it was possible to obtain the Flo- 
quet exponents of different cluster states analytically 
with this technique [21|. By these means it has been 
demonstrated that the unidirectional ring configuration 
of Stuart-Landau oscillators exhibits in-phase synchrony, 
splay states, and clustering depending on the choice of 
the control parameter /?. For /3 = 0, there exists mul- 
tistability of the possible synchronous states in a large 
parameter range. However, when tuning the coupling 
phase to an optimal value 13 = QmT ~ 2Tnn/N according 
to a particular state 777, this synchronous state is monos- 
table for any values of the coupling strength K and the 
time delay r. The main goal of this paper is to find ade- 
quate values of /3 by automatic adaptive adjustment. For 
this purpose, we make use of the speed gradient method 
llSll. which is outlined in the next section. 



Then we evaluate the gradient of uj{x^u^t) with respect 
to input variables: 

V„a;(a;,u,i) = Vu[^ xQ{x,t)Y F{x,u,t). 

Finally, we set up a differential equation for the input 
variables u 



du 
It 



-TVu^{x,u,t), 



(6) 



where F = F^ > is a positive definite gain matrix. 
The algorithm ([6]) is called speed- gradient (SG) algo- 
rithm, since it suggests to change u proportionally to the 
gradient of the speed of changing Q. 

The idea of this algorithm is the following. The term 
— V„aj(x, 77, i) points to the direction in which the value 
of Q decreases with the highest speed. Therefore, if one 
forces the control signal to "follow" this direction, the 
value of Q will decrease and finally be negative. When 
Q < 0, then Q will decrease and, eventually, tend to zero. 

We shall now apply the speed-gradient method to net- 
works of Stuart-Landau oscillators. Since the coupling 
phase /? is the crucial parameter that determines sta- 
bility of the possible in-phase, cluster, and splay states, 
we use this control parameter as the input variable u. 
Setting u = P and x = (zi, . . . , zjv), Eq. ([1]) takes 
the form of Eq.® with state vector x G C^ and in- 
put variable /3 G M, and nonlinear function F{x, P, t) = 
[/(zi), . . . , f{zM)] + Ke^P[Ax{t - r) - x{t)]. 

The SG control equation © for the input variable /3 
then becomes 

f = ~r|.(.,/3,0^-F(|)%..Q(.,0, (7) 
where F > is now a scalar. 



III. SPEED-GRADIENT METHOD 



IV. IN-PHASE SYNCHRONIZATION 



In this section, we briefly review an adaptive control 
scheme called speed- gradient (SG) method. Consider a 
general nonlinear dynamical system 



X = F{x, u,t) 



(3) 



with state vector x € C", input (control) variables u € 
C™, and nonlinear function F. Define a control goal 



lim Q{x{t),t) =0, 



(4) 



where Q{x,t) > is a smooth scalar goal function. 

In order to design a control algorithm, the scalar func- 
tion Q = uj(x,u,t) is calculated, that is, the speed 
(rate) at which Q{x{t),t) is changing along trajectories 
ofEq. @: 

Loix, u, t) = ^%1^ + [V,Q(x-, t)rFix, u, t). (5) 
at 



To apply the SG method for the selection of in-phase 
synchronization we need to find an appropriate goal func- 
tion Q. It should satisfy the following conditions: the 
goal function must be zero for an in-phase synchronous 
state and larger than zero for other states. Hence, a 
simple goal function can be introduced by taking the dis- 
tance of all oscillator phases to a reference oscillator's 
phase Lpi: 



i{x{t),t) 



1 ^ 



^l)^ 



(8) 



fc=2 



Taking the gradient of the derivative along the trajecto- 
ries of the system ([T]) with local dynamics ([2]) one can 
derive an adaptive law of the following form by straight- 
forward calculation. Using Ld{x,P,t) ~ Qi, Eq. ([7]) be- 
comes 
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where we used the abbreviations r^^T = r„(i — t) and 
Vn,T ^ fnit ^ t) for notational convenience. 

Figure [T] presents the results of a numerical simula- 
tion for a random network with N = 6 nodes and unity 
row sum. Throughout this paper wc use F = 1. Ac- 
cording to the numerical simulations decreasing F will 
yield a decrease of the speed of convergence. On the 
other hand, if F is too big, undesirable oscillations ap- 
pear. The model parameters are chosen as in [2l|. In 
Fig. [TJa) it can be seen that the absolute values \zj\ of 
all nodes converge after about 60 time units. Fig. [TJb) 
shows that the phase differences of the different oscilla- 
tors approach zero, which corresponds to the in-phase 
synchronous state. Fig. [Ijc) depicts the evolution of /?. 
The blue dashed line represents the value of the coupling 
phase /3 = Hqt = 0.487r, for which stability was shown 
analytically in [2l[. It can be seen that the adaptively 
adjusted phase comes close to this value. In other words, 
even without knowing the exact values of the system pa- 
rameters, the SG algorithm yields an adequate value of 
P that stabilizes the target state of in-phase synchroniza- 



tion. Fig. [Ijd) shows that the goal function (|S]) indeed 
approaches zero. 

Note that the above choice of the goal function Q is not 
the only possibility to generate a stable in-phase solution. 
Let us consider a function based on the order parameter 



Ri 
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(10) 



It is obvious that i?i = 1 if and only if the state is in- 
phase synchronized. For other cases we have i?i < 1. 
Using this observation wc can introduce the following 
goal function 



N N 



(11) 



j=i 



k=l 



From $ = —T-M^Q2 we derive an alternative adaptive 
law: 
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k=l 3 = 1 
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(12) 



Fig. [2] shows the results of a numerical simulation. As 
before, the amplitude and phase approach appropriate 
values that lead to in-phase synchronization. This time, 
however, the obtained value of /3 does not converge to 
the one for which the analytical approach [l3| has estab- 
lished stability of the in-phase oscillation (blue dashed 
line), but to another limit value. This can be explained 
as follows: There exists a whole interval of acceptable 
values of /3 around the value of the coupling phase for 
which an analytical treatment is possible, such that for 
any value from this interval an in-phase state is stable. 
Our SG algorithm finds one of them, depending upon 
initial conditions. 



V. SPLAY AND CLUSTER STATES 
STABILIZATION 

In this section wc will consider unidircctionally coupled 
rings with N = 6 nodes. That is, the coupling matrix has 
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the following form: 



A 



Let 1 < m < N - I. Then d = LCM(TO,iV)/TO, where 
LCM denotes the least common multiple, is the number 
of different clusters of a synchronized solution. A splay 
state corresponds to rf = A^ while cluster states yield 
d < N. Using similar arguments as those leading to 
Eq. ([5]) we could choose a goal function of the following 
form: 
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N 



J = l 



E ^J ~ '^J^ 



2tt 



(13) 



with j = j mod N. 
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Figure 1: (Color online) Adaptive control of in-phase oscilla- 
tions with goal function Eq. ([8]). (a): absolute values rj = \zj\ 
for j = 1,...,6; (b): phase differences A(j>j = fj — fj+i for 
j = 1, ..., 5; (c) temporal evolution of /?, blue dashed line: ref- 
erence value for f2o = 0.92; (d): goal function. Parameters: 
A = 0.1, a; = 1,7 = 0, ii: = 0.08, r = 0.527r, A^ = 6. Initial 
conditions for rj and ipj are chosen randomly from [0, 4] and 
[0, 27r], respectively. The initial condition for j3 is zero. 
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Figure 2: (Color online) Adaptive control of in-phase oscil- 
lations with goal function Eq. pi|) . (a): absolute values 
^i = l-^j|i (b): phase differences A<j}j — ifij — ifij+i; (c): tem- 
poral evolution of /3, blue dashed line: reference value for 
Qo — 0.92; (d): goal function. Other parameters as in Fig. [T] 



The goal function Eq. P^ has a crucial disadvantage: 
we need to define an ordering of the system nodes. Since 
this is inconvenient for practical applications, we will ex- 



tend the alternative goal function Eq. pT|) such that we 
can stabilize splay and cluster states. First of all, note 
that the following condition holds for splay and cluster 
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state and for other p-cluster states where p are divisors 
of d. 

For example, suppose that the system has six nodes. 
Then states for which conditions P^ and ((TS|) with 
i?c; = 1 for rf = 6 hold are schematically depicted in 
Fig. [3ja),(b),(c). In order to distinguish between these 
three cases, let us consider the functions 



Figure 3: Schematic diagrams of splay (d = 6), three-cluster 
(d = 3), and two-cluster (d = 2) states in panels (a), (b), and 
(c), respectively (A'' = 6). Each cluster contains the same 
number of nodes. 
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Indeed, if we have only three nodes and take Q = 
^ -^^ 6*"^^ Sfc=i e~"^'= as a goal function, we will ensure 
stability of a splay state, as we have verified by numerical 
simulations. Note that this goal function does not need 
a fixed ordering of the nodes. Renumbering all nodes in 
a random way will yield the same goal function. One can 
define a generalized order parameter 
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(15) 



with d e N. However, if we derive a goal function from 
this order parameter in an analogous way as in Eq. dill) , 
this function will not have a unique minimum at the d- 
cluster state because R^ = 1 holds also for the in-phase 
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A splay state (Fig. EJa)) yields /i = /2 = /a = 0, 
while in the 3-cluster state displayed in Fig.|3l^b) we have 
/i = /2 = 0, /a = 1, and in the 2-cluster-state shown 
in Fig. [S^c) /i = /s = 0, /2 = 1. Hence, we obtain 
Sn /p = if and only if there is a state with d clusters, 
where the sum is taken over all divisors of d. x,t) 

Combining all previous results we adopt the following 
goal function: 



04 = 1-/4^) + — E /p(^)' (17) 

p\d,l<p<d 

where p\d means that p is a factor of d. This goal func- 
tion contains fd as the primary contribution for the d- 
cluster state, but also a sum of penalty terms that coun- 
teract reaching other cluster states in which fd is also 
unity. Whenever one of those unwanted cluster states 
is approached, the penalty term will lead to a gradient 
away from it. The prefactor N'^/2 is chosen for con- 
venience to secure faster convergence of the algorithm. 
From /3 = "^'^Qi one can derive the adaptation law 
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(18) 



In Fig. 0] we show the results of a numerical simulation 
for splay state stabilization {d = N = 6, m = 1). The 
phase differences are Acfij = ipj — ^Pj+i = 27r — 2tt/N , 
which corresponds to the splay state. In Fig. UJ^c) one 
can see that the adaptively obtained value of j3 converges 
to that for which stability was shown analytically in |2l| 
(dashed blue line). 

Figures [5] and [6] depict the results of numerical simula- 
tions for two clusters (d = 2, to = 3) and three clusters 
(d = 3, TO = 2,4), respectively. Again we note that the 
obtained value of /3 comes close to the one for which sta- 
bility was shown analytically in j2l| . 

The above results indicate that the speed-gradient 
method is able to drive the network dynamics into the 



desired cluster or splay state by adaptively adjusting the 
coupling phase, where the goal function is chosen accord- 
ing to the corresponding target state. We have, however, 
used only exemplary values of the coupling parameters 
K and t so far. 

For the example of a splay state (4-cluster) in a net- 
work of 4 Stuart-Landau oscillators coupled in a unidirec- 
tional ring we have conducted a more exhaustive analysis 
of the (A', t) plane. Figure [7] shows results in depen- 
dence on the coupling strength K and the coupling delay 
T. According to Ref. [21| there exists an optimal value 
of the coupling phase that enables stability of this state 
for arbitrary values of K and r. We ran simulations 
with 20 different initial conditions chosen randomly from 
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Figure 4: (Color online) Adaptive control of splay state with 
goal function Eq. (|17|) . (a): absolute values Vj = \zj\; (b): 
phase differences A(j>j = ifij — i^j+i; (c): temporal evolution 
of /3, blue dashed line: reference value for ili = 0.96; (d): goal 
function. Other parameters as in Fig. [l] 



Figure 5: (Color online) Adaptive control of 2-cluster state 
(m — 3) with goal function Eq. (|17|l . (a): absolute values rj = 
\zj\; (b): phase differences Acfij — tpj — ipj+i; (c): temporal 
evolution of/?, blue dashed line: reference value for SI3 — 1.08; 
(d): goal function. Other parameters as in Fig. [l] 



the complex interval [—1,1] x [—i,i] for each oscillator 
Zj . Figure [7] shows the fraction fc of those realizations 
that asymptotically approach a splay state after applying 
the speed-gradient method. Wc observe that the speed- 
gradient method is able to control the splay state in a 



wide parameter range. The range of possible coupling 
strengths K does, however, shrink considerably with in- 
creasing time delay r. We conjecture several reasons for 
this shrinking. Firstly, multistability of different splay 
and cluster states is more likely for larger values of K 
and T. which narrows down the basin of attraction for 




(a) 



2jc 



^-.-. 



8ic/6 
4 



(b) 




-0,2 



n^T-lnlZ 



on. 



-0,6 

-0,8 

-1 

■1,2 

-^•\ 
6 

5 

4 

a" 3 

2 





(c) 






\ 11 




^ 


\11 / 




\ / 




\y-y 













(d) 




























^--^_^ 





Figure 6: (Color online) Adaptive control of 3-cluster state 
{m = 2,4) with goal function Eq. (|17p . (a): absolute values 
Tj = l^jl; (b): phase differences A(j}j = ipj — ipj+i, blue dashed 
line: reference value for ^2 ~ 1.03; (c): temporal evolution of 
/?; (d): goal function. Other parameters as in Fig. [T] 




Figure 7: (Color online) Success of the speed-gradient method 
in dependence on the coupling parameters K and t for the 
splay state in a unidirectionally coupled ring of A'' = 4 Stuart- 
Landau oscillators. Other parameters as in Fig. [T] The color 
code shows the fraction of successful realizations. 



VI. CONCLUSION 

We have proposed a novel adaptive method for the con- 
trol of synchrony on oscillator networks, which combines 
time-delayed coupling with the speed gradient method of 
control theory. Choosing an appropriate goal function, 
a desired state of generalized synchrony can be selected 
by the self-adaptive automatic adjustment of a control 
parameter, i.e., the coupling phase. This goal function, 
which is based on a generalization of the Kuramoto or- 
der parameter, vanishes for the desired state, e.g., in- 
phase, splay, or cluster states, irrespectively of the or- 
dering of the nodes. By numerical simulations we have 
shown that those different states can be stabilized, and 
the coupling phase converges to an optimum value. We 
have elaborated on the robustness of the control scheme 
by investigating the success rates of the algorithm in de- 
pendence on the coupling parameters, i.e., the coupling 
strength and the time delay. In this work, we focused on 
the adaptive adjustment of the coupling phase while the 
other coupling parameters were fixed. The input variable 
u in Eq. ([3]) may in general contain all of the coupling 
parameters. Thus, as a generalization, our method might 
be applied to all coupling parameters including the cou- 
pling amplitude and the time delay. In this way control of 
cluster and splay synchronization might be possible with- 
out any a priori knowledge of the coupling parameters. 
Given the paradigmatic nature of the Stuart-Landau os- 
cillator as a generic model, we expect broad applicability, 
for instance to synchronization of networks in medicine, 
chemistry or mechanical engineering. The mean-field na- 
ture of our goal function makes our approach accessible 
even for very large networks independently of the partic- 
ular topology. 
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